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Abstract 

The ,5 = 1 Heisenberg Antiferromagnet is studied in the presence of two kinds 
of local impurities. First, a perturbed antiferromagnetic bond with J' ^ J 
at the center of an even-length open chain is considered. Using the density 
matrix renormalization group method we find that, for sufficiently strong or 
weak J', a bound state is localized at the impurity site, giving rise to an energy 
level in the Haldane gap. The energy of the bound state is in agreement with 
perturbative results, based on s = 1/2 chain-end excitations, both in the weak 
and strong coupling limit. In a region around the uniform limit, J' = J, no 
states are found with energy below the Haldane gap. Secondly, as = 1/2 
impurity at the center of an otherwise even-length open chain is considered. 
The coupling to the s = 1/2 impurity is varied. Bound states in the Haldane 
gap are found only for sufficiently weak (antiferromagnetic) coupling. For a 
s = 1/2 impurity coupled with a strong (antiferromagnetic) bond, no states 
are found in the Haldane. Our results are in good qualitative agreement with 



recent experiments on doped NENP and Y2BaNi05. 
75.10.-b, 75.10.Jm, 75.40.Mg 
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I. INTRODUCTION 



In a seminal paper, Haldane showed that integer spin chains display a gap in the spec- 
trum, where as half-integer spin chains are gaplessS. This was confirmed both by exper- 
imental and numerical calculation^. Currently, the best estimates of the gap for the 
isotropic s = 1 chain are A = 0.41050(2)J, using the density matrix renormalization group 
(DMRG) methodi, and A = 0.41049(2) J using exact diagonalizationi. The dynamical struc- 
ture factor has also been determined in great detail, by numerical studies of finite systems^. 
Most of the experimental effort has in recent years focused on Ni(C2H8N2)2N02(C104) 
(NENP) which show clear evidence of the Haldane gapHi. NENP has not been observed to 
undergo a three-dimensional ordering at any accessible temperature as opposed to CsNiCls 
which yielded the first experimental confirmation of the Haldane gapi. NENP has a fairly 
large single ion anisotropy estimated iohe D / J = O.isi where as CsNiCla is almost isotropic. 
Recently, experiments have been performed on Y2BaNiO This charge-transfer insula- 

tor has long one-dimensional chains of Ni atoms at the center of compressed corner sharing 
oxygen octahedra. Thus, along the chains, the Ni^"*" ions, with s = 1, are separated by 
nonmagnetic O^^ ions. The antiferromagnetic (AF) superexchange interaction between the 
Ni^"*" ions is then mediated through the oxygen sites. Y2BaNi05 has a fairly large exchange 
coupling of J/A;fe=28 5Klli322K0, estimated from high temperature susceptibility measure- 
ments. From inelastic neutron scattering experiments on powdersEl it is known that the 
inter chain coupling is small, J' / J ^ 10^^. Thus, the compound is highly one-dimensional 
and no three-dimensional ordering has been seen. In powder averaged inelastic neutron scat- 
tering experimentslli two gaps have been observed; For fluctuations parallel to the chain, 
A|| = 186 ± 127^^, and for fluctuations perpendicular to the chain, A^ = 99 ± QK. This is in 
good agreement with single crystal measurements of A^ = 100 ± The two gaps are 

explained in terms of a single ion anisotropy D / J = 0.16 ± 0.020, a value close to what is 
found for NENP. 

The one-dimensional Haldane gap systems are by now fairly well understood. However, 
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the important effect of impurities on these systems have only recently been considered. 
Y2BaNi05 offers several possibilities for doping and is thus well suited for such studies. 
Non-magnetic Zn^+ ions can be substituted for Ni^"*", thus effectively severing the Ni chains. 
Experiments have been performed on Zn doped samples, Y2BaNii_2;Zna;O50'0Hl3, confirm- 
ing the picture of chain cutting, with a possibility for an increase in the gap with dopinglii 
due to the fact that the average length of the chain segments is small enough for finite-size 
effects to be observable. Severing of the chains is predicted to have rather dramatic effects 
on the low energy excitations. It can be shown that the AF Heisenberg model is fairly well 
approximated by a valence bond solid (VBS) modelB. The breaking of the valence bonds 
by the introduction of impurities should produce free s = 1/2 spins at the sites next to the 
impurity. The effect of such free s = 1/2 at the chain-ends has been observed in electron- 
spin-resonance (ESR) experiments on NENP by doping with Cu§, thereby introducing a 
s = 1/2 at the impurity site, and by doping with non- magnetic Zn, Cd and Hgi^. These 
findings have recently been contrasted by specific heat measurementJi^ on Y2BaNii_3,Zna,05 
which show a Schottky anomaly consistent with free s = 1. 

Another possibility for doping exists in Y2BaNi05. Carriers, in the form of holes, can 
be added to the system by replacing the off-chain Y^+ with Ca^~^. DiTusa et. al0 show, 
from polarized x-ray absorption spectra, that the added holes are localized on the oxygen 
between the Ni atoms. The addition of a hole on the oxygen sites could effectively turn 
these into s = 1/2 spins, if strongly localized. Ca doping should then have significant 
effects on the spectrum. Powder averaged inelastic neutron scattering have been performed 
on Y2_^.Ca^BaNi05lli showing that Ca doping produces states in the Haldane gap. From 
analysis of the spectral weight per impurity, DiTusa et al.lll have proposed that the magnetic 
disturbance caused by the Ca doping introduces either a triplet or a quartet of subgap bound 
states. 

In order to address the experimental results on doped s = 1 spin chains we therefore 
consider the simple model Hamiltonian 
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L-1 

H=J {^^■^^+l} + ■J E {S^■Si+^} + H,^p. (1.1) 

i=l i=L/2+l 

Here the single ion anisotropy is neglected for simplicity, and S represents s = 1 spins. The 
impurity Hamiltonian ifimp is taken to be either a perturbed bond 

'Sl/2 ■ Sl/2+1, (1.2) 

or a s = 1/2 impurity, in which case we have 

-f^imp = J'^L/2 ■ S' + J'S' ■ Si/2+1, (1-3) 

where S' is the s = 1/2 impurity spin. Schematically the two kinds of impurities are shown 
in Fig. |l]a, |I]b, respectively. For the carrier doped compounds Eq. ( p..l|) may not be a good 
description and a Hubbard like model may be more appropriate. However, in themselves 
the Hamiltonians, Eq. ( |1.1| ), deserves attention and they should form a good starting point 
also for the carrier doped compounds if the carriers are strongly localized. 

Our approach is to study the low-lying excitations in these systems by comparing large- 
scale numerical results, obtained using the density matrix renormalization group (DMRG) 
method, to perturbative results, valid at strong or large coupling, based on effective s = 
1/2 chain-end excitations generated by broken valence bonds. In section || we discuss the 



numerical approach. Section HI addresses the effects of a single perturbed bond at the center 



of a chain, while we in Section |^ present our results for a s = 1/2 impurity. The bulk of 
our results concerns antiferromagnetically coupled impurities, however, in Sec. we briefly 
discuss scenarios for ferromagnetically coupled impurities. 

II. NUMERICAL METHODS 

We use the density matrix renormalization group (DMRG) method proposed by 
whit(liFi. For a detailed discussion of the DMRG method we refer the reader to Refs. |22|J23| . 
In all our runs we truncate each half-space to m = 81 states, and we work in a subspace 
defined by the total z-component of the magnetization, S^, and by the parity with respect 



to reflection around the central bond, or central spin, for even and odd length chains, re- 
spectively. For the DMRG method we therefore denote the different states by S^^. Note 
that, this is the z-component of the spin and not the total spin of the state. 

The idea behind the DMRG method revolves around the notion of a subsystem, in this 
case half of the chain, that is strongly coupled to the universe, here the complete chain. As 
usual in quantum mechanics, the subsystem is then best described by its density matrix, 
whose eigenvalues give the probability for the subsystem to be in a given eigenstate (of the 
subsystem density matrix), under the constraint that the universe (the complete chain) is 
in a given pure state. Here we mean, by pure state, an eigenstate of the Hamiltonian for 
the complete system. It is now clear that if we start to effectively decouple the subsystem 
from the universe, by changing the relevant coupling at the center of the chain, the standard 
DMRG method, where the system is split in half at the impurity site, will become progres- 
sively inaccurate. This can be seen in the limit where the coupling to the universe, i. e. the 
central bond in the chain, is zero, in which case the eigenstates of the density matrix and the 
Hamiltonian for the subsystem are the same and the DMRG method reduces to standard 
real space renormalization which is not comparable to the strongly coupled DMRG in preci- 
sion. We therefore expect to have problems whenever the perturbations that we consider at 
the center of the chain behave in a way that effectively breaks the system in two. However, 
as we shall see, this still leaves a workable region where we can obtain useful results. We 
assume that these difficulties could be overcome by considering more complicated DMRG 
schemes where the system is not split at the impurity site. However, most other schemes 
will reduce the symmetry of the problem and we therefore prefer to work with the standard 
method. 

III. BOND DOPING 

In the case where the Ni based s = 1 Heisenberg antiferromagnets are doped by impurities 
that form divalent non-magnetic ions, such as Zn, Cd, and Hg, we expect the impurities to 
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effectively cut the Ni chains. In the valence bond picturella this leads to free s = 1/2 at the 
neighboring sites. In this section we consider the more general case where the central bond 
in an open chain is different from the rest of the chain, although still antiferromagnetic. The 
Hamiltonian describing the system is Eq. ( |1.1| ) and ( p..2[ ). For the experiments performed 
with non-magnetic impurities the relevant situation is that of an extremely weak bond. This 
can be seen in the experiments of Glarum et al.0 who interpret their results in terms of free 
s = 1/2 chain-end spins. 



A. Weak Coupling 

If the bond at the center of the chain is weak enough we can view the problem as 
two weakly coupled open chains of half the length. At the end of an open s = 1 chain 
we expect to find an effective s = 1/2 excitation, this we shall refer to as a chain-end 
spin. Thus, for sufficiently long chains, the low-lying excitations can be described by an 
effective Hamiltonian of two weakly coupled chain-end spins. This should be valid once 
the coupling of the chain-end spins across the impurity bond is larger than the coupling 
along the chain segment. The latter is known to decrease exponential with the length of 
the chain seg ment§iii, and should thus quickly become negligible. Experimentally this is 
the relevant situation for NENP, since it is difficult to obtain high Zn doping levels^, and 
the average chain segment is therefore likely to be long. Y2BaNi05 offer the possibility of 
higher doping levels. If we denote by the projection onto the s = 1/2 subspace, we 

can write = aS"s=i/2- We need to determine a. Fortunately, an estimate 

can easily be obtained by looking at the expectation value one of the spins at the end of a 
free chain. If we by 1 1 > denote the state that has both the chain-end spins polarized in the 
"up" -state we find 

< i|s^=i|i > ~ < i|p(^/2)s:=iP(^/2)|i > 

= a<iis:j;/2ii>= -. (3.1) 
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Previous workilii have estabhshed that < l|Sf^^|l >^ 0.5320 at the end of a free chain, 
and thus a ~ 1.0640. The effective Hamiltonian can then be written as 

ifeff = aV'S"i/2 ■ S\/2+i J' < J, (3.2) 

where S" denotes a s = 1/2 spin. We then immediately see that the low energy spectrum 
will be given by a ground-state singlet and a higher lying triplet, the splitting between the 
two given by 

= aV, J' < J. (3.3) 

This corresponds to a bound state in the Haldane gap, A. For the isotropic antiferromagnetic 
Heisenberg model that we are studying here we have A = 0.4105J in the absence of any 
impurities. In Fig. ^ we show the energy gap between the two lowest lying states for a 
chain that has a weak bond of J' = J/10 in the middle of an otherwise isotropic chain. 
Since we have two essentially free s = 1/2 chain-end spins at each end of the open chain, 
the ground-state becomes four-fold degenerate. For finite chain lengths we expect these two 
chain-end spins to combine into a ground-state singlet, O"*", with an exponentially low-lying 
triplet, l~iiii. In the thermodynamic limit the triplet and the singlet become degenerate. 
Thus, if we want to calculate the spectrum caused by the impurity we need to look at higher 
excited states. The resulting gap, which remain finite in the thermodynamic limit, can be 
calculated between the 2"*" and 1~ states, where we for computational convenience take the 
1^ state to represent the ground-state since they become degenerate in the thermodynamic 
limit. As is clearly evident in Fig. 0, the gap converges quickly to a value of AE = 0.1141. 
With J' = O.IJ we find that Eq. (|3.3| ) gives AE = 0.1132 in excellent agreement. The good 
agreement between the perturbative results and the DMRG results lends further support to 
the picture of effective s = 1/2 chain-end spins. 

The DMRG method is convenient for studying correlations in real-space. In Fig. |^ 
we show the expectation value of the z-component of the spin along a 100 site chain for 
J' = O.IJ. In the 2+ state we have effectively polarized all four s = 1/2 chain-end spins 
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at sites 1, 50, 51, and 100 in the "up" direction. In order to focus attention on the effects 
of the impurity we have subtracted off < Sf >i-. In the 1~ state the chain-end spins at 
site 50, 51 are combined into a singlet, by subtracting off < S- >i- we therefore cancel out 
the effect of the chain-end spins at site 1 and 100. One sees that the two chain-end spins 
surrounding the perturbed bond have effectively formed a localized state. 

Using the DMRG method we can now map out the energy of this bound state as a 
function of the weak coupling, J'. Due to the problems mentioned in Section |l| it is not 
possible to treat values of J' below J/10 within the standard DMRG. However, since the 
impurity state is localized, the convergence to the thermodynamic limit is very fast as can be 
seen in Fig. |^. Our results are shown in Fig. ^. The filled squares are the numerical results, 
the solid line the Haldane gap for the isotropic chain. A/ J = 0.41050, and the dashed line 
indicates the result from Eq. ( p.3|) . From Fig. ^ it is seen that only for J' < 0.7 J does a 
bound state below the isotropic Haldane gap exist. Around J' = 0.7 J the bound state seen 
in Fig. ^ delocalizes and spreads out over the entire chain. This is shown in Fig. ^ where we 
plot < >2+ — < SI as a function of chain index i for a 100 site chain with J' = 0.7 J. 
Clearly this state is not localized and it resembles what one finds for the isotropic chain@. 
We therefore conclude that it is necessary to weaken the central bond below a critical value 
J'c weak ~ 0.7J in order to create a bound state. 

B. Strong Coupling 

We now turn to a discussion of the opposite limit where the central bond is very strong, 
J' ^ J. For the antiferromagnetic coupling that we consider here, the two neighboring spins 
will be contracted into a singlet. In the limit J' — oo the chain will be severed as was the 
case for a weak bond. We now want to derive an effective Hamiltonian valid for J' ^ J. As 
a starting point it is convenient to consider the case where we only have four spins described 
by the Hamiltonian: 

H = JSi-S2 + J'S2 ■ S3 + JS3 ■ S4. (3.4) 
9 



We do perturbation theory in J/ J' . Note that when J = there is a degneracy since the 
spins, Si and S4 can be in arbitrary states. This is hfted in second order perturbation 
theory in J. The result of second order perturbation theory can be expressed in terms of an 
effective Hamiltonian acting on the low-energy subspace, in which S2 and S3 are projected 
onto the singlet state: 

i/efr = PH'—^—-H'P. (3.5) 

Here 

Ho = J'S2 • S3, (3.6) 
H' = JSi • S2 + JS3 • S4, (3.7) 

and P projects onto the singlet ground-state of Hq, with energy Eq — —2 J. Acting with H' 
on the ground-state of Hq produces the eigenstate of Hq with spin 1 and energy, Ei — — J. 
Thus: 

= ^^(Si • S2 + S3 • S,fP 
2J 

= constant - ^^(Si • S2)(S3 • S4)P 
2 T 

= constant - y < ^ISiS^slQ > S^Sl (3.8) 
where |0 > denotes the spin-0 ground-state of Hq. Using: 

< 0\S^S^^\0 >= < OIS2 • S3IO >= (3.9) 

we obtain: 

^eff = ^^Si-S4. (3.10) 

The same calculation applies to a chain of arbitrary length, i^eff now contains all the 
additional couplings which don't involve the pair of strongly coupled spins, Sl/2 and SLf2+i 
ie.: 
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L/2-2 L-1 ^ j2 

Hes = J X! Sj ■ Sj+i + J ^ Sj ■ Sj+i + -— Si/2_i ■ 5'i/2+2- (3-11) 

i=l i=L/2+2 

This is exactly the same as the original Hamiltonian of Eq. ( p. . 1|) with two sites removed 
from the ends of the chain and J' replaced by 4J^/3J'. Assuming this impurity coupling is 
much smaller than J we write an effective low-energy Hamiltonian by proceeding along the 
lines of Section P^II A] . By introducing spin-1/2 variables, ^^'^ ^"1/2+21 obtain: 

4j2 ^ 

HeS = 3 ^L/2-1 ■ ^L/2+2- (3-12) 

From this effective Hamiltonian we can then derive a simple expression for the bound state: 

4 P 

AE=-a^ — , J':^J. (3.13) 

We have calculated the gap AE = E2+ — Ei- as a function of the coupling, J/ J', to the 
s = 1/2 impurity spin. Our results, extrapolated to the thermodynamic limit, are shown in 
Fig. 1^. The solid line indicates the Haldane gap. A/ J = 0.41050, for the pure chain. The 
solid squares are the numerical results obtained using the DMRG method. The dashed line 
denotes the perturbative result Eq. ( p. 131 ). The open circles indicate a direct diagonalization 



of Eq. ( |3.4| ) with Si and S4 spin one-half. As explained in Section y we are limited to study 
J/ J' > 0.25. For J/ J' = 0.5 the limiting value of the gap is very close to the value of the 
isotropic chain. A, and for J/ J' < 0.5 there appears to be no bound state in the Haldane gap. 
For smaller J/ J' a bound state emerges, and we see good agreement between the numerical 
results and the effective Hamiltonian Eq. (|3.12|) . As was the case for a weak bond we find a 



critical coupling, J/J'^ strong; tiut now with an approximate value of J/J'^ strong ~ 0.5 above 
which no bound state occurs. 

The picture of a quenched singlet in the middle of the chain weakly coupled to two chain- 
end spins on either side, can be substantiated by calculating the expectation value of the 
z-component of the spin along the chain. This is shown in Fig. ^ for L = 100, J/ J' = 0.25. 
Again we have subtracted off the 1~ state configuration. The quenched singlet is clearly 
visible at the center of the chain and we see a large signal from the two polarized chain-end 
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spins at i = L/2 — 1 and i = L/2 + 2. Since one site has effectively been removed from each 
of the two half-chains Fig. ^ appears inverted with respect to Fig. ^. At the point where 
the bound state disappear we again find that the expectation value < >2+ — < >i- 
spreads out over the entire chain. 

We may think of a magnon as behaving, in many respects, like an ordinary particle of 
mass m = A/v'^, where v is the spin- wave velocity. Thus it may, at first, seem surprising 
that no bound state exists for a weak perturbation of the coupling, J' J since there is a 
well-known theorem of one-dimensional quantum mechanics which states that an attractive 
potential always produces a bound state, no matter how weak it i s§. Our result is not 
surprising if one looks in a bit more detail at the mapping of the spin problem into an 
effective one-particle quantum mechanics problem. 

We represent the spin operators at long length-scales in terms of the non-linear a (NLcr) 
model field, (f) (obeying the constraint = 1) and the associated canonical momentum, 
1 = {l/vg)(p X d(p/dt, as: 

S, ^ s(-l)'(/) + l. (3.14) 
The Hamiltonian density becomes: 

Acting on the vacuum, cf) creates the triplet of single-magnon states. A local perturbation 
in the exchange couplings leads to additional local terms in the field theory Hamiltonian 
density of the form: 

5n = vVix)(^^^ +K'(x)l2, (3.16) 

where V{x) and V{x') are some short-range functions. Note that we do not get simply 
V{x)4>'^, since this term is trivial due to the constraint, 0^ = 1. 

We sometimes approximate the NLcr model by a free massive field theory model with 
Hamiltonian densityil: 
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7^^^^^ + ^f^V + -«^^ (3-17) 

2 2\dxJ 2v^ ' ^ ' 

where the constraint on <p is now relaxed and 11 is the associated cononical momentum. 
This gives the dispersion relation: 



E = ^A^ + {vky. (3.18) 

(We have shifted the wave- vector by tt.) We may easily treat the additional impurity term, 
vV {x){d(f) / dxY . This gives the Schrodinger equation: 

- (1/2X - [V{x)^Ij']' = mEtl). (3.19) 

We obtain this equation for each of the three magnon polarizations. Note that this does not 
correspond to a standard potential scattering term. We could obtain a potential scattering 
term from a term in the NLcr model Hamiltonian of the form V{x){(l)^Y or V{x)l^. However, 
these terms break the SU (2) symmetry and therefore can't appear in the Heisenberg model. 
The latter term corresponds to a local magnetic field. The former could arise if we have a 
crystal field term in addition to the Heisenberg exchange. 

The peculiar derivative potential term does not necessarily lead to a bound state if it is 
weak. To demonstrate this we consider the case of a step function: 

V{x) = (|a;| > a) 

V(x) = -^0 {\x\ < a). (3.20) 

Solutions to the Schrodinger equation will have discontinuous ip' at x = ±a. Defining ip^a"^) 
to be the limit of ip'{x) as x — > a from above or below respectively and similarly for %l)'{—a^) 
we obtain: 

v.V) = V'V)[i-2H] 

V^'(-a-) = V''(-a+)[l - 2V^]. (3.21) 

The important thing to notice is that if Vq < 1/2, the derivative does not change sign at 
the discontinuity. There can then be no solution which vanishes at x ^ ±oo. This follows 
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since, assuming an even wave-function, it must have the form cosh[y 2|i?|m/(l — 2Vo)x for 
|x| < a which is increasing as we move away from 0. It has to turn into exp[— Y/2|ii^|m|x|] 
for |x| > a. This can't happen without a change in sign of the derivative, [ip itself must 
be continuous.] Alternatively, we could assume an odd wave-function but we have the same 
problem. If Vq > 1/2, there is a solution. It is an odd function of x; a sine function for 
|x| < a. The binding energy is determined by: 



2\E\m 



tan^^^^a=^2K,-l. (3.22) 

This has a solution for all Vq > 1/2. As Vq — 1/2+, we get: 

\E\m ^ {2Vo - lY/2a\ (3.23) 

This provides a counter-example to the statement that there is always a bound state for 
a weak derivative potential term. On the other hand, there is a rigorous proof of the fact 
(see, for example Ref. ^) for an ordinary potential term. This argument is admittedly quite 
incomplete since we haven't analyzed the other, \^(a;)F term, nor dealt carefully with the 
passage from the NLcr model to the free model. (This additional term seems to give a type 
of non-linear Schrodinger equation.) Ignoring these complications, and assuming that the 
step function result is universal, the above argument suggests that as J' J'^ the critical 
coupling where the bound state disappears, the binding energy behaves as: 

Ef, oc (J' - J'f. (3.24) 

This is not inconsistent with our numerical results, but substantially more numerical work 
for J' close to would be needed in order to verify this dependence. 

Our findings in this section are in disagreement with a recent Schwinger-Boson calcula- 
tionii. 

IV. 5 = 1/2 IMPURITIES 

We now turn to a discussion of s = 1/2 impurities in a s = 1 AF Heisenberg chain. 
Experimentally two scenarios have been studied. First, Cu^+ (s = 1/2) was substituted for 
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Ni^"*" (s = 1) in NENP0. Here, the introduced s = 1/2 impurity structurally takes the same 
place as the original Ni ion. Experimentally one finds a weak ferromagnetic coupling to the 
neighboring chain-end spins^. Secondly, carriers were added to Y2BaNi05 by substituting 
Ca for Y. DiTusa et al.0 argue that the holes become localized on the O sites between the 
Ni ions. Here the situation is different since the introduced s = 1/2 impurity does not take 
the place of a s = 1, instead, it is placed between two Ni ions. Thus, the superexchange 
between the two Ni ions is broken and presumably replaced by a direct exchange to the 
s = 1/2, J', which is likely to be large compared to J. Aharony et ahS have argued that 
for La2_a;(Sr,Ba)2;Cu04_j, the analogous coupling should be large, \J'\ ^ \J\. However, 
the appropriate sign and value of J' for a simple linear chain Heisenberg model capable of 
describing the experiments on Y2_a;Ca3;BaNiO50 is an open question. In the following we 
shall limit our discussion to the case of an antiferromagnetic coupling, J' > 0. 



A. Weak Coupling 

Following the approach in Section [111 A|, for a weak bond at the center of the chain. 



we can try to find an effective Hamiltonian for a weakly coupled s = 1/2 impurity at the 
center of the chain. In this case we again view the sites L/2 and L/2 + 1 of the s = 1 chain 
as being occupied by chain-end spins that effectively behave as if they were spin one-half. 
By projecting onto the s = 1/2 subspace we represent the chain-end spins at these two 
sites by aS", where S" is now a real spin one-half. We then obtain an effective low energy 
Hamiltonian. With S' representing the s = 1/2 impurity spin we have: 

= aJ'[S"L/2 ■ S' + S' • S"i/2+i]. (4.1) 

The spectrum for this effective Hamiltonian consists of a ground-state doublet and an excited 
doublet and quadruplet. The ground-state doublet has 5*^ = 1/2"*", and energy E = —aJ'. 
The excited doublet, = 1/2^, has energy E = 0. Finally, the quadruplet, = 3/2^, 
has energy aJ' /2. From this effective Hamiltonian we see that we should have two low-lying 
levels with energies in the Haldane gap with energies given by 
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aJ , 



(4.2) 



and 



AE3/2+ = -aJ'. 



(4.3) 



For a weak ferromagnetic coupling we would have obtained A£^i/2 
3aJV2. 



aJ'/2, AE1/2+ 



In order to observe these two levels with the DMRG method we need to determine the 
relevant quantum numbers for the complete system. We begin by considering the case where 
the two half-chains surrounding the s = 1/2 impurity are of even length. We consider the 
situation where the chain-end spins at sites L/2 and L/2 + 1 are coupled weakly with a 
coupling g, g J', to the chain-end spins at the sites 1 and L, respectively. Since the chain 
segments have an even length this coupling, g, is antiferromagneti^. It is now straight 
forward to obtain the level diagram shown in Fig. ^a. If the two half-chains have an odd 
length the coupling, g, is ferromagneti^ and we obtain the level diagram shown in Fig. §b. 
Note that, since we in either case have an even number of valence bonds we do not get a 
change of parity from the rest of the chain. Thus, we see that in the thermodynamic limit 
the ground-state is eight-fold degenerate with an eight-fold and sixteen-fold degenerate level 
above it. In our implementation of the DMRG method we cannot use the total spin as a 
good quantum number but only the total z-component. In order to determine the two levels 
given by Eqs. ([4.2|) and ( [4.3|) we then see from Fig. || that we need to calculate the energy 
differences -E3/2- —E-^/2+ and -E5/2+ — -^'3/2+- Only for L/2 odd is 3/2+ the true ground-state. 
The Haldane gap to bulk excitations (i. e. excitations not localized at the center and/or 
ends of the chain) will be between the 5/2^ and 3/2+ states. 

In Fig. ^ we show the energy gap, AE{L), as a function of chain length L, for a s = 1 
chain with a s = 1/2 impurity. The coupling to the s = 1/2 impurity is J' = 0.2 J. The 
symbols denote the gaps: E^i2-{L) — E-i/2+{L) (triangles) to the state 5/2~ (Haldane gap), 
E^/2+{L) — E^i2+{L) (squares) to the state 5/2+, and E^i2-{L) — £'3/2+ (i^) (circles) to the 
state 3/2"^. Clearly, we have two low- lying excitations well separated from the Haldane gap. 
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For the 3/2" state we find the gap to be 0.2073 J, compared to 0.2128 J from Eq. 
Although at large L we have minor problems with the numerical results for the 5/2+ state 
we can also in this case estimate the gap to ~ 0.31 which compares favorably to the result 
0.3192 from Eq. ( pl) . 

It is also possible to study the spatial extension of the bound states. However, as opposed 
to Section |T|, < S'^ > for the state 3/2+ is non-zero across the impurity site. It is therefore 
difficult to give a quantitative interpretation of the results since the subtraction of < >3/2+ 
will produce more complicated effects than just the cancellation of the contribution from the 
free chain-end spins at the sites 1 and L. However, one can still differ between extended and 
localized states. We found that the two states 3/2" and 5/2+ clearly were localized around 
the impurity site for J' = 0.2 J whereas the state 5/2" was extended. 

We can now calculate the low- lying excitations as a function of coupling to the s = 1/2 
impurity. Our results are shown in Fig. [1^. The symbols denote the gaps: -E5/2- — Es/2+ 
(triangles) to the state 5/2" (Haldane gap), -E5/2+ — (squares) to the state 5/2+, 

and -E'3/2- — -E3/2+ (filled squares) to the state 3/2". The sohd line is the Haldane gap, 
AE = 0.41050 J. The dashed line is the perturbative result Eq. (f4.2|), and the dot-dashed 
line Eq. ( |4.3| ). For small J' we see excellent agreement with the perturbative results. For 
J' < J'^ ^ 0.45J a bound state, calculated from the 3/2" level, appears, and for J' < 
J'c 2^ 0.35J also bound state, calculated from the 5/2+ level, is stable. 

B. Strong Coupling 

As before, we try to arrive at an effective Hamiltonian capable of describing the low 
energy excitations. Again we shall assume that the coupling to the s = 1/2 impurity spin, 
J', is always positive corresponding to an antiferromagnetic coupling. As a starting point 
we consider the three central spins 8^/2,8^/2+1 and 8', the s = 1/2 impurity. Disregarding 
the rest of the chain, these three spins will be described by a Hamiltonian 

H,^p = J'[Sl/2-S' + S'-Sl/2+i] (4.4) 
17 



The spectrum of -ffimp is as follows: The ground-state is a quadruplet = 3/2"'" with 
E = -3/2 J', then a doublet = 1/2- with E = -IJ', another doublet = 1/2+ with 
E = 0, a. quadruplet = 3/2" with E = J'/ 2 and finally a sextuplet = 5/2"'" with 
energy E = J'. Since the ground-state of -ffimp is a multiplet with s = 3/2, the effective 
Hamiltonian becomes that of an effective s = 3/2 impurity in the otherwise uniform AF 
s = 1 chain. If we by P^^/^) denote the projection onto the s = 3/2 ground-state subspace 
of i^imp, we can write 

p(3/2)g^^2p(3/2) ^ /5S;'^3/2. The effective Hamiltonian can then be 

written: 

L/2-2 



L-1 

i=l i=L/2+2 



(4.5) 



where the term describing the s = 3/2 impurity is given by: 

=3/2 "T ^s=3/2 " '-'L/2+2J- (4.6) 



PJ[Sl/2-1 ■ ^=3/2 + K=3/2 ■ Sl/2+2] 



Here 8^^3/2 denotes the effective s = 3/2 spin. We need to determine p. We begin by writing 
down the matrices for Sl^2 Sl^2+i projected onto the 3/2"'" ground-state subspace of 
iJimp. These have to be proportional to a spin-3/2 representation. We find: 



< 3 1 ^2 I > < "^is*^ 1"^ > 

^1 L/2\ c2^\ L/2-\-l\ 
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(4.7) 



Clearly, this is just (3/5)5*^^3/2, as expected. Thus, if we project 8^/2 and 8^/2+1 onto the 
3/2"'" subspace we get in both cases a result proportional to a spin-3/2 by a factor 3/5, i.e. 
(3/5)8^=3/2. Thus, 



/3 = 3/5. 



(4.8) 



The coupling to the s = 3/2 impurity is therefore fixed at J3/5. 
Let us consider the more general case where we instead of Eq. 



have 



H,^^{s = 3/2) = K[Sl/2-i ■ S'U/2 + S'U/2 ■ Sl/2+2], (4.9) 

i. e. we regard the coupling to the s = 3/2 impurity as a free parameter, K. We would 
expect bound states to appear for K < For very small K, K <^ Kf'/'^\ we could 

again assume that effective s = 1/2 chain-end spins on the sites L/2 — 1 and L/2 + 2 would 
be a valid description, and we would arrive at an effective low-energy Hamiltonian consisting 
only of the two chain-end spins, and the s = 3/2 impurity: 

ifefr = o^K[Sl„_, ■ S;U/2 + S;U/2 ■ S^/2+2], K « (4.10) 

Here, a ~ 1.064, and ^'1/2+2 ^ ~ spins. The low-energy Hamiltonian, 

Eq. ([4.10|) , has the following spectrum. A ground-state doublet with energy E = —aK5/2, 
an excited quadruplet at E = —aK, another quadruplet at = 0, and finally a sextuplet 
aX E = aK3/2. For sufficiently small K all three bound states should be stable giving rise 
to energy levels in the Haldane gap given by the corresponding energy differences, aK3/2, 
aKh/2, and aKA. 

Returning to the full effective Hamiltonian, Eq. ( [4.5| ), the question of whether bound 
states appear comes down to the question of whether J3/5 is below or above the threshold, 
Kf'/'^\ As we shall see below, the numerical results indicate that J3/5 is too large for any 
bound states to be stable. It is then not possible to write down an effective low-energy 



Hamiltonian based on the s = 1/2 chain-end spins, as we did in Eq. (|4.10|) , and it is 
necessary to work with the full effective Hamiltonian, Eq. ([4.5|) , which, due to the size of 
the corresponding Hilbert space, is not very useful. 

We have calculated the gap £'5/2+ — -£^3/2+ for J/ J' = 1, 0.2, and 0.1 using the DMRG 
method. For large J' this is the lowest-lying excitation. The results are shown in Fig. |T^. In 
this case the convergence to the thermodynamic limit is rather slow and we have difficulty 
extrapolating the results to a high accuracy, we therefore show the unextrapolated results 
for a 120 site chain. Here the solid squares are the DMRG results, and the solid line denotes 
the isotropic Haldane gap of A = 0.41050J. In this case the DMRG method can be used 



for even very large J'. As is clearly evident from Fig. [Tl| we never observe any states in 
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the Haldane gap. A calculation of £"3/2- — £"3/2+ also didn't show any bound states in the 
Haldane gap. We therefore conclude that the coupling to the effective s = 3/2 impurity, 
J3/5, is too large for any bound states to appear. 

We have also calculated the expectation value of the z-component of the spin for the 
values of J / J' = 1, 0.2, and 0.1. Our results are shown in Fig. |T2l The circles are the 



results for J/ J' = 1, the squares for J/ J' = 0.2 and the triangles for J/ J' = 0.1. In order to 
differ between extended and localized states we have subtracted < Sj >3/2+- As opposed to 
Section |T| the spin configuration that we subtract, here < Sf >3/2+, is non-zero across the 



impurity making it difficult to give a quantitative interpretation of these results. However, 
clearly the state is in all cases extended over the whole length of the chain. 

Thus, we conclude that for a strong antiferromagnetic coupling to the s = 1/2 impurity 
no bound states exists below the Haldane gap. The fact that DiTusa et ahlll observes such 
states, indicates that, the appropriate (AF) J', for a simple linear chain Heisenberg model 
description of their results, should be rather small. 

V. FERROMAGNETIC COUPLINGS 

We have only performed the numerical DMRG calculations for antiferromagnetically 
coupled impurities. It would be interesting to repeat the calculation using ferromagnetic 
couplings. Due to the rather large amount of CPU time required we have not done this. 
Since we in our implementation of the DMRG method do not work with the total spin as 
a good quantum number, but only with it's z-component, ferromagnetic interactions are 
generally also more difficult to treat, due to the fact that the ground-state belongs to a 
different multiplet for even and odd L/2. However, the effective Hamiltonians we have 
obtained in the previous sections can be derived for a ferromagnetic coupled impurity in a 
straightforward way. Below we briefly discuss what one finds. 

For a central bond that is weak, \J'\ ^ |J|, and ferromagnetic in an otherwise anti- 
ferromagnetically coupled s = 1 chain, we expect the results of Sec. [Ill A| to still apply. 
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In particular Eq. ( p.3|) should still be correct. However, in this case the ground state of 
Eq. (|]2D is the triplet. 

For a central bond that is strong, \J'\ ^ \J\, and ferromagnetic, the two central spins, 
Sl/2 and Sl/2+i, will have a s = 2 ground-state. We then obtain an effective Hamiltonian 
in terms of an effective s = 2 impurity: 

L/2-2 L-1 
i/eff = J E S.-S,+i + J E S,-S,+i + iJi^p(s = 2), (5.1) 

i=l i=L/2+2 

where the term describing the s = 2 impurity is given by: 

H,mp{s = 2) = J[Sl/2-i ■ V^^^Sl/2V^^^ + V^^^Sl/2+iV^^^ ■ SL/2+2] 

= lJ[^L/2-l ■ K=2 + K=2 ■ SL/2+2] ■ (5.2) 

Here 7 is determined by the projection of 8^/2 and 8^/2+1 onto the s = 2 subspace, 
V^'^'^Sl/2V^'^^ = 78'/^2- It is easy to see that 7 = 1/2, thus, the effective coupling to the s = 2 
impurity is J/2. If the coupling to the s = 2 impurity is regarded as a free parameter we 
would expect bound states to appear below a threshold, Kj,^\ However, in the absence of 
any numerical data, we have no way of telling whether J/2 is above or below such a threshold 
and we cannot draw any conclusions. Let us briefly consider the case where the coupling to 
the s = 2 impurity is assumed to be a free variable K. For K K^^ we can again use the 
picture of effective s = 1/2 chain-end spins to write an effective low-energy Hamiltonian: 

Hes = «-^[S'^/2-i ■ ^"=2 + ^"=2 ■ ^"l/2+2\- previous sections, 8'^/2-i 8'^/2+2 '^'^^ 

effective s = 1/2 spins. This low-energy effective Hamiltonian has the following spectrum: 
A ground-state triplet with energy E = —aK3, followed by a quintuplet at E = —aK, 
another quintuplet at = 0, and finally a septuplet a.t E = aK2. For sufficiently small 
K, K <ti Kf \ we then expect bound states in the Haldane gap to occur at energies, 2Ka, 
?)Ka, and bKa. The effective Hamiltonian, Eq. ( [5.1| ), corresponds to = J/2, since this 
is not a small coupling, a low-energy Hamiltonian based on the s = 1/2 chain-end spins is 
presumably not valid. Since we do not know K'f'\ it is not possible to tell if bound states 
occur. 
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If a s = 1/2 impurity is coupled with a weak ferromagnetic bond to the otherwise 
antiferromagnetically coupled s = 1 chain, we expect bound states to occur. As already 
mentioned in Sec. [IV A| , we expect these to have energies AE1/2- = aJ'/2, AE1/2+ = 
aJ'3/2, with a ~ 1.064. 

For as = 1/2 impurity coupled with a strong ferromagnetic bond to an otherwise uniform 
s = 1 chain, we know from Sec. [IV B| that, the ground state of -ffimp, Eq- (|4^ ) should be a 
sextuplet, = 5/2^^. If we by P^^/^) denote projection onto the 5/2 subspace, we can define 
a constant, 6, by V^^^'^^Sl/2V^^^'^^ = 6S"^^/2- We can then write an effective Hamiltonian in 
terms of an effective s = 5/2 impurity in the following way: 

L-l 

ifefr = J E S.-S,+i + J E S,-S,+i + ifi^p(s = 5/2), (5.3) 

«=1 i=L/2+2 

where the term describing the s = 5/2 impurity is given by: 

H,^^{s = 5/2) = J[Sl/2-i ■ V^'/''>Sl/2V'^'/'^ + V^'/^^Sl/2+iV^"''^ ■ SL/2+2] 

= 5J[Sl/2-1 ■ ^'1=5/2 + ^'1=5/2 ■ Sl/2+2]- (5.4) 

In this case S = 2/5. The effective coupling to the s = 5/2 impurity is then 2 J/5. If we 
regard this coupling as a free parameter, K, we expect bound states to occur for K below a 
threshold, i^f 

are not able to determine if 2 J/5 is below such a threshold. As above, 
let us consider the case where the coupling to the 3 = 5/2 impurity is a free variable K. 
For K < JsTf /2) we can again use the picture of effective s = 1/2 chain-end spins to write 
an effective low-energy Hamiltonian: iJefr = ' ^"=5/2 + ^"=5/2 ' '^l/2+2]- above, 

'^L/2-1 ^^'^ ^L/2+2 effective s = 1/2 spins. This low-energy effective Hamiltonian has 

the following spectrum: A ground-state quadruplet with energy E = —aK7/2, followed by a 
sextuplet ai E = —aK, another sextuplet ai E = 0, and finally an octuplet at E = aK5/2. 
For sufficeintly small K we then expect bound states in the Haldane gap to occur at energies, 
aK5/2, aK7/2, and aK6. The effective Hamiltonian, Eq. ( [5.3[ ), corresponds to i^' = 2J/5, 
this coupling cannot be regarded as small and we are not in the regime K < Thus, 
the description based on the s = 1/2 chain-end spins is of limited use. However, whether 
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2 J/5 is below or above K^^l'^^ is an open question. 

We see that only in specific cases is it straight forward to derive a useful low- energy eSec- 
tive Hamiltonian. For a strong ferromagnetic bond, and for a s = 1/2 impurity coupled with 
a strong ferromagnetic bond, further numerical work will be needed in order to determine 
whether bound states in the Haldane gap exists. 

VI. DISCUSSION 

We have shown that, by the introduction of an AF bond with J' ^ J, J' > in an 
otherwise isotropic AF s = 1 chain, a triplet bound state below the Haldane gap occurs for 
sufficiently large and small J'. The energy of the bound state is in very good agreement 
with perturbative results based on effective s = 1/2 chain-end excitations created by broken 
valence bonds. In a region 0.7 J ^ J'^ ^^^-^ ^ J' ^ J'c strong ~ 2.0J no bound state is found 
to within numerical precision. If a s = 1/2 impurity is introduced with an AF coupling 
J' ^ J, J' > in an isotropic AF s = 1 chain we find that, only for J' below i ~ 0.45J do 
bound states occur. For J' < J'^ ^ ^ 0.45J a bound state appears, and for J' < 2 ~ 0.35J 
a second bound state is stable. At large impurity couplings, J' > J, we find no bound states. 

In experiments on Y2_:cCaa;BaNiO50, states below the Haldane gap have been found. 
In a simple linear chain spin model of Y2_x.Ca3;BaNi05, our results would imply that the 
coupling to a s = 1/2 impurity, introduced by the Ca doping, should be rather small, if 
it is antiferromagnetic. This may seem surprising since we would expect \J'\ > |J| to be 
the more likely scenario in Y2_xCaxBaNi05. However, the coupling could be large and 
ferromagnetic in which case we cannot rule out the existence of bound states. Secondly, 
the non-zero single ion anisotropy term, D{S^y, relevant for Y2_2,.Caa;BaNi05 could create 
bound states in the gap even for a large antiferromagnetic coupling to the s = 1/2 impurity. 
Also, the substitution of Y^+ by Ca^+ will lower the crystal field symmetry on close by 
Ni sites, justifying new local anisotropy terms in a spin Hamiltonian. Possibly such terms 
could create bound states. DiTusa et al.0 argue that the magnetic disturbance caused by 
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the Ca doping introduces either a triplet or a quartet of subgap bound states. Our results 
indicate that several different multiplets of bound states are likely to be present, some of 
which have spin larger than 1/2, even for the simple scenario of Ca doping giving rise to 
s — 1/2 impurities, thus justifying the large contribution to the spectral weight from the 
impurities seen in the experiments. 

Our finding that a critical coupling is needed to form bound states in this problem 
may seem surprising at first, since one would expect that even an infinitely small impurity 
potential would form a bound state in one dimension. However, in a NLcr model description 
one can argue that the introduction of impurities do not give rise to simple local potentials 
due to the constraint |</)p = 1. This might change in a model where the rotational symmetry 
is broken by the introduction of a single ion anisotropy term, D{Sf)^. Such a term is relevant 
for the description of NENP and Y2BaNi05. It would therefore be very interesting to repeat 
the above DMRG calculations with D y^O. 
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FIGURES 

FIG. 1. Schematic picture of the two different perturbations. In a) is shown the perturbed 
bond with couphng strength J' at the center of a s = 1 chain with couphng J. In b) is shown the 
s = 1/2 impurity coupled with couphng strength J' to the rest of the s = 1 chain. 

FIG. 2. The energy gap AE{L) = E2+ (L) — Ei- (L) as a function of chain length, L, for a weak 
bond, J'/ J = 0.1. The limiting value of the gap is 0.11405. 

FIG. 3. The bound state in real space for a chain with L = 100. < Sf >2+ — < Sf >i- for a 
weak bond, J' / J = 0.1. Clearly visible is the localized state around the impurity site. 

FIG. 4. The energy gap, /S.E = E2+ — Ei- (solid squares), as a function of the weak perturbed 
bond, J' /J. The solid line is the isotropic Haldane gap, A/J = 0.41050. The dashed line indicates 
the perturbative result, Eq. (|3.3D . 

FIG. 5. < 5f >2+ - < 5f > I- as function of chain index i for a weak bond, J'/ J = 0.7, close 
to the critical coupling. 

FIG. 6. The energy gap AE = E2+ —Ei~ as a function of the strong perturbed bond, J/ J' . The 
solid squares are the DMRG results. The solid line is the isotropic Haldane gap, A/J = 0.41050. 
The dashed line indicates the perturbative result, Eq. ( p.l3| ). The open circles denote exact diag- 
onalization results of Eq. (|3.4D . 

FIG. 7. The bound state in real space for a chain with L = 100. < Sf >2+ — < Sf >i- as a 
function i for a strong bond, J/ J' = 0.25. The quenched singlet right at the center of the chain is 
clearly visible. 

FIG. 8. The low- lying levels for a weakly antiferromagnetically coupled s = 1/2 in a spin-one 
chain. 



28 



FIG. 9. The energy gap, AE{L), as a function of chain length L, for a s = 1 chain with a 
s = 1/2 impurity. The couphng to the s = 1/2 is J' = 0.2 J. The symbols denote the gaps: 
E^I2-{L) — £'3/2+ (L) (triangles) to the state 5/2" (Haldane gap), E^i2+{L) — £'3/2+ (i) (squares) 
to the state 5/2+, and E^i2-{L) — E^i2+{L) (circles) to the state 3/2~. 

FIG. 10. The low-lying excitations as a function of coupling, J/ J', to a weakly coupled s = 1/2 
impurity. The symbols denote the gaps: £5/2- —-£'3/2+ (triangles) to the state 5/2" (Haldane gap), 
~ (squares) to the state 5/2"^, and £3/2- — -£'3/2+ (filled squares) to the state 3/2~. 

The filled squares denote the DMRG results for a 120 site chain. The solid line is the Haldane gap, 
AE = 0.41050 J. The dashed line is the perturbative result , Eq. ( [4.2]) , and the dot-dashed line 
Eq. (P). 

FIG. 11. The gap, £5/2+ ~-^3/2+) ^ ^ function of coupling, J/J', to the s = 1/2 impurity. The 
filled squares denote the DMRG results for a 120 site chain. No extrapolation to the thermodynamic 
limit has been performed. The solid line is the Haldane gap, AE = 0.41050J. 

FIG. 12. < Sf >5/2+ — < >i/2+ as a function i, for J/J' = 1 (circles), J/ J' = 0.2 (squares), 
J /J' = 0.1 (triangles). 
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Fig 2. Sorensen, Affleck 
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Fig 3. Sorensen, Affleck 
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Fig. 4. Sorensen, Affleck 
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Fig. 5. Sorensen, Affleck 




Fig. 6. Sorensen, Affleck 
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Fig. 7. Sorensen, Affleck 
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Fig. 8. Sorensen, Affleck 
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Fig. 9. Sorensen, Affleck 
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Fig. 10. Sorensen, Affleck 
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Fig. 1 1 . Sorensen, Affleck 
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